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The observed leptonic mixing pattern could be explained by the presence of a discrete flavour 
symmetry broken into residual subgroups at low energies. In this scenario, a residual generalised CP 
symmetry allows the parameters of the PMNS matrix, including Majorana phases, to be predicted in 
terms of a small set of input parameters. In this article, we study the mixing parameter correlations 
arising from the symmetry group A 5 including generalised CP subsequently broken into all of its 
possible residual symmetries. Focusing on those patterns which satisfy present experimental bounds, 
we then provide a detailed analysis of the measurable signatures accessible to the planned reactor, 
superbeam and neutrinoless double beta decay experiments. We also discuss the role which could be 
played by high-precision measurements from longer term projects such as the Neutrino Factory. This 
work provides a concrete example of how the synergies of the upcoming experimental programme 
allow flavour symmetric models to be thoroughly investigated. Indeed, thanks to the rich tapestry of 
observable correlations, we find that each step of the experimental programme can make important 
contributions to the assessment of such flavour-symmetric patterns, and ultimately all patterns that 
we have identified can be excluded, or strong evidence found for their continued relevance. 
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I. INTRODUCTION 

The existence of three families of fermions in the Stan¬ 
dard Model (SM), identical in all properties apart from 
their masses, is as yet unexplained by any physical prin¬ 
ciple or mechanism. Moreover, the discovery that both 
quarks and leptons permit complementary but distinct 
descriptions in terms of the flavour states which cliag- 
onalise the weak interactions and the states which di- 
agonalise their mass terms has shown that the connec¬ 
tion between families betrays a precise structure which 
is an essential component in our description of the phys¬ 
ical world. Explaining the origins of this flavour struc¬ 
ture has been a recurring theme in proposed extensions 
of the SM. One such programme is the application of 
discrete flavour symmetries, predominately in the lepton 
sector, where the flavour quantum numbers are associ¬ 
ated with a new symmetry and particles are assigned to 
its irreducible representations. This can provide a way 
to unify the three families into a single mathematical ob¬ 
ject. However, as the lepton masses are known to be 
distinct, any non-abelian symmetry can only be exact 
above the scale of mass generation. Nevertheless, its ex¬ 
istence at high energy shapes the theory, and the residual 
symmetries which survive the breaking procedure at low 
energies can play an important role in the structure of 
flavour observables. 

The paradigm of a non-abelian flavour symmetry 
breaking into residual symmetries has been used by many 


* peter.ballett@durham.ac.uk 
t silvia.pascoli@durham.ac.uk 
t jessica.turner@durham.ac.uk 


authors to make predictions about the six mixing an¬ 
gles and phases which constitute the Pontecorvo-Maki- 
Nakagawa-Sakata (PMNS) matrix and parameterize neu¬ 
trino mixing: $ 12 , $ 13 , O 23 , the Dirac phase <5 and two 
Majorana phases 021 and 0 : 31 - For a recent review of 
such models see e.g. Ref. HQ. Although many of the ear¬ 
liest models were designed to predict a very small value of 
013 , a prediction now firmly ruled out GH, a number of 
models remain consistent with the current data, many of 
which are based on groups taken from the A( 6 n 2 ) family: 
A(96) 13El, A(150) H, A(600) [5] and A(1536)[5]. This 
connection was strengthened in Ref. [10] which showed 
that, based on only a few generic model building assump¬ 
tions, if the full PMNS matrix is to be specified by the 
symmetry alone (so-called “direct” models CD. the only 
possible predictions which would agree with current data 
are those arising (minimally) from A( 6 n 2 ). In general, 
however, the existence of residual symmetries amongst 
the leptonic mass terms may not fully specify the mixing 
pattern. In these “semi-direct” models [I] the symme¬ 
tries reduce the degrees of freedom necessary to describe 
the mixing parameters by defining a correlation between 
previously independent parameters (sometimes known as 
mixing sum rules HIHll)- Often these correlations be¬ 
tween mixing angles and phases can be derived from quite 
generic analyses of the residual symmetries present in a 
system without needing to specify a full UV-complete 
theory [T5U2T1 . As such, focusing on these relations can 
be an effective way to compare a wide class of models to 
data HZH3]. 

To date, the fine structure of the PMNS matrix has 
been inaccessible to experiment, preventing the study of 
subtle parameter correlations. Current measurements of 
the mixing angles have 3er uncertainties of around 6.9% 
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on 0 i 2 , 7.3% on 0 13 and 16.7% on 0 23 EEjj All three 
of the CP phases are unconstrained at this significance 
level; although, some low-significance hints for a maxi¬ 
mally CP-violating value of the Dirac phase 6 ~ 3tt/2 
have been observed [25, (2Bj. Therefore the upcoming 
experimental work will focus on two key topics: the pre¬ 
cision determination of the mixing angles and the first 
measurements of the CP phase S. This will open the door 
for studies of the correlations between mixing angles, and 
between mixing angles and phases, that are predicted by 
models of flavour symmetries. 

In the most popular formulation of models with dis¬ 
crete flavour symmetries, the constraints on the mass 
matrices used to derive the PMNS matrix cannot remove 
a number of complex phase degrees of freedom. This re¬ 
sults in an inability to predict the Majorana phases and, 
in general, lessens the predictivity of the model. How¬ 
ever, by imposing a generalised CP symmetry (GCP), 
phase information may be accessible and dictated by the 
flavour structure itself. This can lead to very predictive 
scenarios, where all 6 mixing parameters are related to a 
small number of input parameters [27]. GCPs were first 
explored in the context of discrete and continuous groups 
in Ref. |28j ; however, they have recently been revived due 
to the question of consistency between a CP symmetry 
and discrete flavour group [271 [21 [M]. This has lead to 
interesting work studying the predictions of models with 
imposed flavour and GCP symmetries for a number of 
groups such as A 4 [21 [31] , S 4 [ZHE2H35], A(48)[2S[37] 
and A (96) [35] along with more comprehensive analy¬ 
ses of the families A(3n 2 ) and A( 6 n 2 ) [52MT] . (See also 
Ref. [42] and Ref. (43 for further applications of GCP 
symmetries.) 

In this article, we present a detailed analysis of a sin¬ 
gle group: the alternating group on 5 elements, A 5 . This 
was first introduced to leptonic flavour physics in Ref. [44j 
via the study of Golden Ratio mixing 21 SS] : a possible 
pattern of the PMNS matrix with $13 = 0, 0 2 3 = 7 t /4 
and a value of 0 i 2 related to the golden ratio p = 1+ 2 V ^ , 
tan 0i 2 = 1/ip. This pattern has been shown to be a 
prediction of a number of different models based on A 5 
[4411451 - H1 : however, it is not the only fully-specified mix¬ 
ing pattern associated with direct models based on this 
group. If a Z 3 subgroup is preserved among the charged 
leptons and a Z 2 x Z 2 is preserved among the neutrinos, a 
pattern with 0 i 3 = 0, 0 2 3 = 7 r /4 and cos0i 2 = ip/\[3 can 
be found EDEBEI]. Further patterns are also possible 
when a Z 2 x Z 2 symmetry is preserved in the charged- 
leptons whilst a different Z 2 x Z 2 remains in the neutrino 
sector, which predict a large value of 0 i 3 , 0 i 3 ~ 17.9° 

El- 

Needless to say, the patterns above are in severe ten¬ 
sion with the current global data by dint of their 0 i 3 
predictions alone; however, the possibility remains that 


symmetries such as these do not completely survive at 
low energies and that a semi-direct approach may remain 
viable. In this work, we consider the group A 5 with a 
GCP symmetry, deriving the most general GCP transfor¬ 
mation which could be implemented for this group. We 
assume that the flavour group with GCP is broken into a 
set of residual symmetries at low energies insufficient to 
fix all of the oscillation parameters. We compute all pos¬ 
sible predictions for the induced correlations amongst the 
mixing parameters. These are compared to the current 
data, and we identify those which are compatible with 
the current bounds. The viable patterns that we iden¬ 
tify are highly predictive, expressing all six parameters of 
the PMNS matrix in terms of a single unphysical angle. 
We take particular care in assessing the phenomenology 
of the correlations between mixing angles and phases for 
these viable models: discussing their accessibility to reac¬ 
tor, long-baseline and neutrinoless double beta decay ex¬ 
periments, and highlighting particularly interesting sig¬ 
natures to be tested. Although we have restricted our 
attention to the group A 5 , this can be seen as an illus¬ 
trative choice and we would like to stress the rich but 
moreover readily testable phenomenology which exists in 
the residual symmetry framework, much of which arises 
from the predictions taken as a whole instead of resting 
on single generic types of parameter correlation. 

The work presented in this paper is divided into two 
main parts. In Section [TT] and Section HI we discuss 
the assumptions behind our framework, and explain the 
technical steps in our derivation of the correlations. In 
Section mi we focus on the predictions themselves, pre¬ 
senting some simplified formulae for the correlations be¬ 
tween observable quantities and identifying their most 
interesting phenomenological signatures. We also study 
a number of ways that the correlations can be tested 
by upcoming reactor, superbeam and neutrinoless dou¬ 
ble beta decay experiments, as well as possible longer- 
term experiments such as neutrino factories. We make 
our concluding remarks in Section [V] 


II. RESIDUAL FLAVOUR AND GENERALISED 
CP SYMMETRIES 

We assume the presence of a discrete flavour symme¬ 
try, G, at high energies. To unify the three flavours, we 
assume that the fields are assigned to a 3-dimensional 
irreducible representation of this group, and the general 
multiplets of leptons 'F transform as 

(K a -A p(g)af}'&ft > 

where p : G —> GL(3, C) represents a unitary representa¬ 
tion of 43 Ensuring the existence of the 3-dimensional 


1 For alternative global analyses of oscillation data, see Ref. m ■ 


2 Throughout this paper we will assume that all representations 
are unitary and therefore all group elements are represented by 
unitary matrices. 
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irreducible representation p restricts us to non-abelian 
groups. 

As neutrinos are known to oscillate, they cannot have 
degenerate masses and therefore the non-abelian flavour 
group, G, cannot be a symmetry of our low-energy ef¬ 
fective lagrangian. Therefore, we assume that the full 
flavour symmetry must be broken at low energies into 
two abelian residual symmetry groups, G e and G„, which 
are unbroken in the charged-lepton and neutrino sectors, 
respectively. We denote the leptonic mass terms in the 
low-energy effective theory by 

G\v — (cl)o (^a)ck/5 (cr)/ 3 d~ (cr)c* (m^)a/3 (cl)/3 

+ ^l) a {m v ) a p{u L ) p + ^L) a (ml) af5 {ul) 0 . 

where the Greek indices are flavour indices and the mass 
matrices are 3x3 and complex-valued. Due to the anti¬ 
commutation of the fermionic fields, one can show that 
the Majorana mass matrix m v is restricted to being com¬ 
plex symmetric. 

We assume that there exist residual symmetries acting 
on the left-handed charged- and neutral-leptons. If we 
denote a general element of these subgroups by g e £ G e 
and £ G„, the fields transform according to the fol¬ 
lowing rules 

(e L )a ->■ P {9e) af 3 (e L )/3, and (u L ) a -* p {g„) a0 {y^)p. 

Combining these relations with the lagrangian in Eq. ([T]) 
leads us to matrix relations which the mass terms must 
satisfy if the residual symmetries are to be preserved at 
low energies 

m x m{ = p(g e y (m x m[)p{g e ) , (2) 

m v = p (g„) T m v p (g„). (3) 

These relations constrain the forms of the mass matrices 
and, as we will show, knowledge of their existence can be 
used to derive a form of the PMNS matrix, ITpmns- 

Working from a bottom-up perspective, we would like 
to deduce the phenomenological consequences of a given 
choice of residual symmetries G e and G„. The choice 
of residual flavour groups is constrained in two ways. 
Firstly, G e and G v must be subgroups of the unbroken 
group G. Secondly, the possible residual flavour symme¬ 
tries must be subgroups of the largest symmetry allowed 
by the mass terms in Eq. ([!]). These maximal symmetries 
are best identified in the basis where both mass terms 
are diagonal. In this basis, the most general symmetry 
of the charged lepton mass matrices is C/(l) 3 , which has 
discrete subgroups of the form G e = Z m for any m or a 
direct product of such groups. For the neutrino residual 
symmetry, the argument changes due the assumed Ma¬ 
jorana nature of the mass term. In this case the largest 
possible symmetry is smaller, Z 2 x Z 2 , which leaves us 
with only two choices, G v = Z 2 or G„ = Z 2 x Z 2 . 

In addition to a flavour symmetry G, we also assume 
the presence of a GCP symmetry. A CP symmetry is 


understood as a combination of charge-conjugation and 
a parity transformation; however, a generalised CP sym¬ 
metry is one which also acts on the flavour indices whilst 
making this transposition. In Ref. [25], it was shown that 
ensuring the consistency of a discrete flavour symmetry 
and a CP symmetry often requires the introduction of a 
non-trivial generalised CP symmetry. 

We define our generalised CP symmetry to act on a set 
of fields ’F a as 

d'a — t A a( 3'I ’p, 

where X a p is assumed to be a unitary matrix so as to 
preserve the kinetic terms in the lagrangian and \k c de¬ 
notes the conventional CP conjugate appropriate for the 
Lorentz representation of the field fk. If a discrete flavour 
symmetry is present, G, then a generalised CP symmetry 
must satisfy a consistency equation EH ESI 

Xp{g)*X * = p(g'), (4) 

where g and g' are elements of G. For a faithful rep¬ 
resentation p, this relation can be seen as establishing a 
mapping from g to g' which preserves the structure of the 
group and therefore defines a group automorphism. In 
Ref. [30], it was pointed out that a physical GCP trans¬ 
formation must be restricted to a single irreducible rep¬ 
resentation, and is related to a class-inverting automor¬ 
phism of G, meaning that g' is mapped to an element 
in the conjugacy class of <? _1 . We restrict our consider¬ 
ation to involutory GCP transformations, requiring that 
the application of the transformation twice is equivalent 
to the identity, and therefore the X matrix satisfies an 
additional constraint 

XX* = 1. (5) 

As with the flavour symmetry G, for our GCP symme¬ 
try to leave the lagrangian invariant the mass matrices 
must satisfy further constraints. The GCP symmetry ex¬ 
changes the hermitian conjugate terms in the lagrangian 
of Eq. <[T|) , which remains invariant if the mass matrices 
obey the relations 

X T m„X = m*, ( 6 ) 

X*(mAin^)X = ( m\m \)*. (7) 

If Eqs. ([ 6 ]) and ([7]) are unbroken relations at low en¬ 
ergies, it can be shown that all CP violating effects of 
the PMNS matrix vanish E3(M]. However, if only one 
of these relations is preserved, the consistency of the 
flavour and CP symmetries leads to novel constraints on 
the PMNS matrix. In this paper, we assume that the 
GCP symmetry is broken in the charged-lepton sector 
but is preserved in the neutrino sector. For this to be 
consistent, the X matrix must map the elements of the 
neutrino residual symmetry to themselves, 

Xp(g v )*X* = p(g v ). 


( 8 ) 
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In summary, we assume a discrete flavour symmetry G 
and a GCP symmetry implemented by X at high energy 
scales which is assumed to break into a subgroup G e act¬ 
ing on the charged-lepton mass terms and another sub¬ 
group G„ which along with the GCP symmetry acts on 
the neutrino mass terms. This leads to a system of con¬ 
straints which the mass matrices must satisfy: Eq. ([2]), 
Eq. (J3|, Eq. (|6j) and Eq. ([8]). In the next section we will 
show how knowledge of these constraints alone can be 
used to predict the PMNS matrix, including its Majo- 
rana phases. 


A. Constructing the PMNS matrix using 
symmetry constraints 

Constraints on the mass matrices, as derived above, 
lead to restrictions on their allowed form, and subse¬ 
quently to the matrices required to diagonalise them. We 
focus first on the charged leptons. A constraint on the 
charged-lepton mass matrix of the form in Eq. (|2|) can be 
rephrased as a statement of commutation 

\p(g e ), ( m x m{)\ = 0. 

As the matrix p(g e ) is unitary and the matrix H = m\m\ 
is hermitian, there exists a basis such that they are si¬ 
multaneously diagonalised, 

3U e s.t. U\U e = 1, p(g e )d = ulp{g e )u e , H d = UlHU e , 

where p(g e ) d and H d denote diagonal forms of the ma¬ 
trices p(g e ) and H. As the charged leptons have distinct 
masses, H is full rank. This implies that U e is unique 
(up to re-phasing and re-ordering of columns). If p{g e ) 
is also known to be full rank, U e can also be found by 
diagonalising this operator. In this way, by insisting on 
the relation in Eq. pi), we can compute U e solely from 
the group element p{g e ), and the symmetry alone speci¬ 
fies the mixing matrix. However, a complication arises if 
p(g e ) is not full rank. In this case, p(g e ) does not have a 
unique diagonalising matrix, as in any basis in which it 
takes diagonal form, further SU(2) transformations can 
be performed freely in its degenerate eigenspace. With¬ 
out knowledge of the mass matrix, our knowledge of p(g e ) 
will only allow the identification of the family of diagonal¬ 
ising matrices of p(g e ), and U e must take a more general 
form 

u e = u 0 R e (4>, 7 )$, 

where Uq is any matrix which diagonalises p(g e ), -R e (<A> j) 
is a complex rotation in the degenerate subspace of p(g e ) 
by an angle <fi with a phase 7 and $ is a diagonal matrix 
of phases. 

In the neutrino sector, we have three constraints to 
consider on the mass terms: one from the flavour sym¬ 
metry, one from the GCP symmetry and one ensuring 


their consistency. Under a change of flavour basis, the 
matrix X is mapped to 

A' -> rfXU*. 

As X is unitary and symmetric, Takagi factorization al¬ 
lows us to express it as X = for some unitary matrix 
f2 which implies that we can choose a basis where A' be¬ 
comes trivial m- In fact, this basis is not unique and 
the remaining freedom can be used to further diagonalise 
p(g u ). In this basis the constraint in Eq. ([6]) implies that 
the mass matrix is real-valued 

€ R. 

As p(g^) is diagonal and commutes with this matrix, we 
know that the mass matrix in this basis must be diagonal 
up to a basis change in the degenerate subspace of p(g v ). 
As it is purely real, the most general additional basis 
transformation required to bring it into diagonal form is 
a rotation in 2-dimensions 

u„ = nR u (e), 

where 9 is the angle describing the real rotation. There 
remains the possibility that the diagonal mass matrix is 
not positive definite, in which case a diagonal re-phasing 
must occur. Without further knowledge of the mass ma¬ 
trix this cannot be predicted, and in consequence, the 
Majorana phases can be predicted only up ±7r or multi¬ 
ples thereof 

We see that the GCP symmetry in the neutrino sector 
has specified a special basis in which the residual flavour 
symmetry elements are diagonal and the mass matrix is 
real. In this way, GCP symmetries help to fix some of the 
re-phasing degrees of freedom associated with diagonal¬ 
ising matrices and allow for the prediction of Majorana 
phases. Combining the results for the charged-lepton sec¬ 
tor and the neutrino sector, we find the full PMNS matrix 
is given by 

U P MNS = ^Re^,l)U i o nR^9), (9) 

where R e and R v denote two unspecified rotations ( R e = 
1 if ord(<? e ) > 2). We make two further simplifications: 
$ is removed by re-phasing the charged leptons, and we 
note that the angles 9 and <j> need only be defined over 
the interval #,</>€ [0,7r), as shifts by 7r can be absorbed 
by unphysical redefinitions of the complex phases. 

III. MIXING PATTERNS FROM A s 

The preceding section showed how the assumed resid¬ 
ual flavour and GCP symmetries can lead to expressions 


! We will always work with the Particle Data Group parameteri¬ 
zation of the PMNS matrix ED, in which the Majorana phases 

( : a 21 : Q 31 \ 

l,e 2 ,e 2 ) which 

take physical values on the intervals a.ij E [0,27r). 
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for the PMNS matrix. In this section, we will derive the 
possible mixing matrices which arise by this method for 
the group A 5 . First we will discuss the structure of A 5 
and the possible subgroups eligible to be taken as resid¬ 
ual symmetries. Then we will derive the form of the 
most general GCP transformation. In the subsequent 
subsections, we consider all viable combinations of CP 
and residual subgroups and present those patterns which 
are consistent with the current global data [25] , 

A. Subgroups and GCP symmetries 

The group A 5 can be defined as the group of even per¬ 
mutations on 5 elements. It has the abstract presentation 

{S, T \ S 2 = T 4 5 = (ST) 3 ), 

where S and T are the generators of the group and all 
group elements can be expressed by a word made from 
these distinguished elements. The structure of this group 
and its representation theory have been discussed in the 
physics literature before (see e.g. Refs. [44] and |48|) and 
we abstain from deriving the explicit representations, de¬ 
ferring the reader to these references instead. However, 
we will briefly recap those features of the group and its 
representations most pertinent to our subsequent analy¬ 
sis. 

We assume that the lepton doublets are assigned to 
a 3-dimensional representation. A 5 has two distinct 3- 
dimensional irreducible representations, and in the fol¬ 
lowing we will always work with the representation 3 of 
Ref. |3Hj. We have checked that the final results do not 
change if we choose the alternative 3-dimensional rep¬ 
resentation instead. For our chosen representation, the 
generators S and T can be expressed by 


- 1 

0 

°\ 

1 / 

' 1 


0 

-1 

0 

and T = - 

v 

-<Pg 

v° 

0 

V 

2 \ 

-Vs 

1 


where = 1 ^' 2 V ^ represents the Galois conjugat^j of 
if. We note that this is a real representation and in our 
chosen basis all group elements are real; it also forms a 
subgroup of SU(3). 

The non-identity elements of A 5 are either of order 2, 3 
or 5 and can be partitioned into 4 conjugacy classes: one 
of order 2 elements (15 members), one of order 3 elements 
(20 elements) and two of order 5 elements (12 members 
each). The centre of A 5 is trivial, and the identity alone 
forms one additional conjugacy class. 

As described in Section [Tl] to compute the PMNS ma¬ 
trix from residual symmetries, we must first find U e , the 
matrix which diagonalises the generator of the residual 


4 In this case Galois conjugation exchanges the two solutions of 

the minimal polynomial over the rationals x 2 — x — 1 = 0. 


symmetry of the charged-lepton mass term. To do so, we 
must identify the eligible residual symmetry groups for 
G e . This symmetry must be an abelian subgroup of A 5 , 
of which there are 4 kinds (up to isomorphism). Three of 
these subgroups are cyclic groups: Z 2 , Z 3 and Z 5 . These 
are the groups generated by a single element g, and its 
members are those powers of g less than or equal to its 
order 

( 9 ) = {g n I s.t. n < ord(g)}. 

As any element can be taken as the generator of a cyclic 
group, there are 15 distinct subgroups of Z 2 in A 5 , 10 of 
Z 3 and 6 of Z 5 . The diagonalizing basis for these groups 
is simply the basis which diagonalises the generator p(g). 

In addition to the cyclic subgroups, there are also non- 
cyclic abelian subgroups in A 5 . These are isomorphic 
to the Klein four group Z 2 x Z 2 which is generated by 
distinguished pairs of order -2 elements, 

{91,92 \gl =g\ = (gig2) 2 )- 

In fact, the 15 order-2 elements in A 5 can be divided 
into 5 triplets which (with the identity) define distinct 
four-element groups. For these non-cyclic groups, the 
diagonalizing basis is defined as that which diagonalizes 
the two generators simultaneously. 

Therefore, the different choices for the residual symme¬ 
try of the charged-lepton mass term can be divided into 
four categories depending on the preserved subgroup: 
G e £ {Z 2 ,Z 3 ,Z 5 ,Z 2 x Z 2 }. For the residual symmetry 
of the neutrino Majorana mass term, we are restricted to 
taking subgroups of Z 2 x Z 2 , leaving us with two options: 
a single Z 2 or the full Klein group Z 2 x Z 2 . 

B. Deriving X 

The matrix X which implements the generalised CP 
symmetry, as discussed in Section[TT] must satisfy XX* = 

1 and be related to a class-inverting automorphism of the 
group, 

Mg £ A 5 , 3h g £ A 5 s.t. (X*p(g)X)* = p(h- 1 )p(g- l )p{h g ), 

where p is our chosen irreducible representation, gener¬ 
ated by the matrices S and T. We shall derive the most 
general form of X for the group A 5 by exploiting our 
knowledge of the automorphism structure of the group. 
The automorphism group of A 5 is S 5 (see e.g. [52]). 
and we identify two important subgroups: inner and 
outer automorphisms. The inner automorphism group, 
Inn(A 5 ), comprises those automorphisms which can be 
represented by conjugation by a group element, 

<kh € Inn(A 5 ) Mg £ A 5 , <j> h (g) = h -1 gh. 

This group can be found by considering the map from el¬ 
ement (h £ A 5 ) to inner automorphism (<fih(g) = h~ 1 gh ), 
and applying the first isomorphism theorem, 

Inn(As) = A.§ / Z (A. 5 ) = A 5 , 
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where the final step uses the fact that A 5 has trivial cen¬ 
tre, Z{A.§) = 1. Therefore, the inner automorphisms 
of A 5 are given by A 5 itself. The outer automorphism 
group is defined as the quotient of the full automorphism 
group by the inner automorphism group. For A 5 it fol¬ 
lows from our discussion above that this is the unique 
group of two elements Z 2 . Our derivation of X is greatly 
simplified by A 5 being an ambivalent group, where each 
element is conjugate to its inverse. For such groups, the 
class-inverting automorphisms are also class-preserving. 
All inner automorphisms of a group are class-preserving, 
but the two properties are not equivalent as there do ex¬ 
ist class-preserving outer automorphisms for some groups 
|53l [54] . However, for the case of A 5 we have a single non¬ 
trivial outer automorphism to check, and this automor¬ 
phism maps elements of order 5 from one conjugacy class 
to the other. Therefore, in the present case, we conclude 
that the class-preserving automorphisms are precisely the 
inner automorphisms. 

We can therefore simplify our defining constraint on 

A, 

3h € A s , Vff € As (X*p(g)X)* = p{h- 1 )p(g)p(h ), 

where the element h is the same for all elements g. As 
we are working with a real representation, we can always 
change basis so that all group elements are given by real 
matrices, and we use this fact with Eq. © to make fur¬ 
ther simplifications 

V ff GA 5 Xp(g)X* =p(h- 1 )p(g)p(h), 
which is equivalent to a commutation relation, 
[p(h)X,p(g)\ = 0. 

We can then invoke Schur’s lemma to infer that as p(h)X 
commutes with all the elements of an irreducible repre¬ 
sentation, it must be a scalar matrix: p(h)X = Al, for 
some complex constant A. Requiring that XX* = 1 con¬ 
strains p(h 2 ) = 1/1 A| 2 . However, by closure the element 
on the left must be a member of A 5 and, as our repre¬ 
sentation is unitary, we conclude that A is just a complex 
phase, A = e 16 for 0ei. Therefore, h must be an order 
2 element, and the most general form of X which imple¬ 
ments an involutory class-inverting automorphism for A 5 
is given by 

A' = e w p(h) s.t. ord (h) = 2. 

In our basis, the consistency relation in Eq. © implies 
that the A' matrix must commute with the generator S of 
the residual Z 2 symmetry in the neutrino sector. There¬ 
fore not all choices of h can be consistently implemented, 
and there will be only 3 non-trivial X matrices (up to 
global phases) for any given S. These are the three ele¬ 
ments of the Klein four group associated with S. If we 
work in the basis where this group is diagonal, we find 
that 

A ° ° \ 

X 1 = e w 0 -1 0 , 

\° 0 -V 


and A 2 and X 3 can be defined as permutations of this 
matrix, where the row of the positive entry is denoted by 
the subscript. It is necessary for us to find the basis in 
which X is trivial. The necessary change of basis is given 
for Xi by f]j, where 

(l 0 0 \ 

Hi = e ie/2 0 i 0 , 

V° 0 V 

and similar definitions hold for i = 2 and * = 3. In this 
basis, the generator of the residual symmetry in the neu¬ 
trino sector is diagonal and the GCP action is trivial. 
However, due to the degenerate subspace in the Z 2 gen¬ 
erator, it remains possible that the neutrino mass matrix 
is only block diagonal and requires an orthogonal trans¬ 
formation to fully diagonalise it. This rotation must be 
in the plane of the degenerate subspace for the matrix S. 
Given these two elements, the most general form of the 
matrix which maps between neutrino flavour and mass 
bases is given by 

u v = nR(9), 

where R(0) is an orthogonal matrix effecting a rotation 
in either the 12-, 13- or 23-plane. We also note at this 
point that the overall phase included in our definition of 
H can be seen to have no physical effect, and will be set 
to zero in what follows. 

In the following subsections, we consider all possible 
residual symmetry groups in the charged-lepton sector. 
In the first three sections we consider the charged-lepton 
residual symmetry to be given by each of the 10 Z 3 sub¬ 
groups, 6 Z 5 subgroups and 5 Z 2 xZ 2 subgroups. For each 
of these subgroups, the basis is found which diagonalises 
its elements, U e . Due to the order of these subgroups, 
this diagonalizing matrix is predicted exactly with no re¬ 
maining degrees of freedom ( e.g. R e = 1 in Eq. ©)■ The 
PMNS matrix is then constructed combining U e with one 
of the forms of U v found above by consideration of the 
residual CP symmetry, 

bfpMNS = Ul^lR{9). 

Finally in Section fill F[ we consider less restrictive sym¬ 
metries when U e is not fully specified by symmetry alone 
{Re / 1 )- 

For each configuration considered in this section, the 
arbitrariness in eigenvector ordering and phasing is ac¬ 
counted for by considering all permutations of rows and 
columns of the PMNS matrix. From these permuted ma¬ 
trices, we compute the mixing angles and phases, and 
these are compared to global data. We report all patterns 
of mixing angles found by this process which are consis¬ 
tent with the current 3<r regions as reported in Ref. [25] . 

C. Predictions from G e = Z3 and G„ = Z 2 

When the residual symmetry in the charged leptons is 
taken as Z 3 , the diagonalising matrix of the residual sym- 
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FIG. 1. Mixing angles for Z 3 as a function of the internal 
parameter 6 . This pattern predicts | sin <5| = 1 and sin cm = 
sinQ 3 i = 0. The shaded regions show the 3<r allowed region 
for the corresponding mixing angle according to current global 
data |25l . 


metry generator T is uniquely specified (up to diagonal 
rephasings and permutations). We have considered the 
10 Z 3 subgroups of A 5 which could act as the residual 
symmetry of the charged leptons. Although many dif¬ 
ferent group elements lead to viable mixing patterns, all 
viable solutions can be described by a single matrix after 
a suitable permutation and redefinition of the unphysical 
parameters. This leads to a single viable set of correla¬ 
tions between the angles. The angles can be derived from 
the PMNS matrix, 


( ^75 0 \ 

^ PMNS= -75 73 75 P 13(0) ’ (10) 

V 76 73 71/ 

where i?i 3 (0) denotes a rotation in the 13-plane by an 
angle 0. This leads to the following expressions for the 
mixing angles 

1 2 
sin 2 0 i 2 = -— . 2 , sin 2 d 13 = - sin 2 0, 

3 — 2 sm 0 o 

sin 2 023 = \ ■ 

This pattern is continuously connected to the tribimaxi- 
mal mixing pattern [55] which is recovered at 0 = 0, and 
is an explicit example of a trimaximal pattern ]5H ] where 
\U a2 \ = l/y/z V« £ {e,/x,r}. We have plotted these 
mixing angle predictions for the full range of the unob¬ 
servable parameter 0 in Fig. |T] in which the coloured re¬ 
gions show the current 3<r global intervals from Ref. :[25| . 

The Dirac phase for this pattern depends discretely on 
the value of 0. It can be shown that 


whilst the Majorana phases can be shown to take CP 
conserving values { 0 : 21 , 0 : 31 } 7 { 0 , 7 r} for all values of 0 . 
From Fig. [l] we see that there are two intervals of the 
unphysical parameter which lead to mixing angles which 
satisfy the current global 3cr bounds. Due to the symme¬ 
try of the expressions, these two solutions offer identical 
predictions for values of 0 3 2 , 0i 3 and 0 23 ; however, one 
of these values lies in a region with <5 = f whilst the 
other predicts 6 = ^7 Therefore, there are two sets of 
predictions from the order-3 elements, differing only in 
their prediction for the Dirac CP phase. 


D. Predictions from G e = Z5 and G v = Z2 

If the residual charged-lepton symmetry is assumed to 
be G e = Z 5 , there are 6 subgroups which we need to 
consider which could lead to a viable set of mixing pa¬ 
rameters. We find two distinct sets of correlations which 
are viable for some range of the unphysical parameter, 
leading to three distinct sets of mixing angle predictions. 

The first set of predictions can be derived from the 
following matrix, 


0 


I y/2 + (£ y/2~\~tp 

Cpmns = - 75 I R 13 ( 8 ), 

\ 1 V 1 

\ \/ 4+2y \/ 4+2c p \/2y 

which leads to mixing angles expressed by 


( 11 ) 


1 


sin 2 0 i 2 = ——j 

1 + cos 2 0 

sin 2 0 23 


. 2 n sm 
sm 0i 3 = — 


9 ' 


2 °23= X- 


The mixing angle predictions from this pattern are shown 
on the left of Fig. [2] as a function of the unphysical pa¬ 
rameter 0. In this case, the Dirac phase is maximally CP 
violating with cos 6 = 0. However, as with the order-3 
elements, the sign of sin 5 depends on the parameter 0, 


6 = 



0 & (0, f), 

0 <£ (f ,tt), 


and the Majorana phases are again given by CP con¬ 
serving values, {a 2 i,a 3 i} 7 { 0 , 7 r}, although the precise 
values cannot be determined in this framework. 

The second viable pattern arising from Z 5 is shown 
on the right panel of Fi g. [2] This can be derived from 
a matrix similar to Eq. (Ill but distinct in the relative 


phasing between the columns, 


Cpmns — 


0 


vT+2 Tp V 4 +2 V 75 
1 y 1 

7 4T 2ip 7 4+2 y 72 y 


775 *13(0)- (12) 


6 = 



0 s (0, |), 

0 <£ (f ,tt), 


This relative phase difference, which arises from the 
choice of alignment between the matrix implementing the 
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FIG. 2. The two patterns of mixing angles for Z 5 as a function of the internal parameter d. Both patterns predict sin 021 = 
sin 031 = 0. The right pattern predicts | sin<5| = 0 whilst the left pattern predicts | sind| = 1. The shaded regions show the So- 
allowed region for the corresponding mixing angle according to current global data 125] . 


GCP symmetry X and the generator S of the residual Z2 
symmetry of the neutrino mass term, crucially affects the 
mixing angle 023 and leads to mixing angles which can 
be expressed by 


sin 2 012 


! 


1 + ip 2 cos 2 0 ’ 


sin 2 0i3 = 


sin 2 0 

TT^ 2 ’ 


sm- V 23 = 2 


^ ^sin 0 + yj\ + ip 2 cos 0^ 


1 + ip 2 cos 2 0 


All CP phases take CP conserving values for this pattern 
of mixing parameters. The precise value of the Dirac 
phase again depends on 0, 


6 = 



9 £ (Oj ?)> 
0 e (f ,tt). 


The appearance of CP conservation can be explained 
as, although it was not imposed explicitly, the gener¬ 
alised CP symmetry remains accidentally unbroken in 
the charged-lepton sector. This second pattern leads to 
two distinct allowed intervals in 0 once we restrict the 
mixing angles to lie in the current 3 cr intervals. Due to 
the symmetry of the curves, these two viable sets of mix¬ 
ing angles are distinguished only by their predictions for 

$23- 

Both of the above mixing patterns are continuous ex¬ 
tensions of the well known GR mixing pattern ( a.k.a. 
GR 1 or GRA) [ 45 } [ 46 ], which is found at 0 = 0 . If we ex¬ 
pand in the small parameter r = \[2 sin 0i3 ffH, we find 
that both patterns arising from Z5 lead to the prediction 
for s = \/3 sin 0 12 — 1 }TT] given by, 


s = \l —- 1 + \l 3 r -+ 0 (r 4 ). ( 13 ) 

V 2 + ip \ 2 + ip 4 ' ' 

The two patterns are distinguished by their predictions 
for 023; however, both can be seen as subcases of a more 


general model predicting the 0i2 correlation in Eq. ( 13 ) 
and further obeying a linearised atmospheric sum rule 
first derived in Ref. m, 


a = ^ r cos <5 + 0 (r 2 ,a 2 ). 

From this relation, it is clear how the maximal angle of 
the first pattern a = 0 is associated with the vanishing of 
cos S , while the contrasting non-trivial predictions of 023 
in the second pattern is due to the CP-conserving value 
of 5 , sin S = 0 . However, the correlation of maximal at¬ 
mospheric mixing and CP conserving values of S is not 
an artefact of linearisation and holds exactly, as is shown 
in the right panel of Fig. [ 2 ] We shall consider these cor¬ 
related maximal predictions as a measurable signature in 
Section HVBI 


E. Predictions from G e = Z2 x Z2 and G v = Z2 

The only non-cyclic Abelian subgroup in A 5 is the 
Klein four group Z 2 x Z 2 . If we take this as the resid¬ 
ual symmetry in the charged lepton sector, we find two 
patterns of mixing angles which differ only by their pre¬ 
dictions for 023 - Consistent predictions exist for all of the 
3 cr range of 0i3, and the predictions for the other mixing 
angles can be seen in Fig. [ 3 ] 

The first pattern can be derived from the following 
mixing matrix, 


1 1 

f ^ 

% 


1 

0 

° \ 

t^PMNS = 77 

<Pg 

1 


° 

COS0 

sin0 

Z 1 


-V 5 

<PgJ 

v° 

— sin0 

cos e) 


From this matrix, we find the mixing angles can be ex- 
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pressed by, 


1 + ip g [cos 2 0 + sin (29)] 

3 — (p g [sin 2 9 — sin (20)] 

1 + ifg [sin 2 9 — sin (20)] 

4 ’ 

1 + ip [cos 2 0 — sin (20)] 

3 — (fig [sin 2 0 — sin (20)] 

For this pattern the complex phases are given by CP 
conserving values: sine) = 0 and { 0 : 21 , 0 : 31 } Q {0,7r}. 
The true value of S can be shown to depend on 0 

5 _ f 0 31.7° < 0 < 58.3° or 121.7° < 0 < 159.1°, 

— ( 7 T else. 


sin 2 0 i 2 

sin 2 0 i 3 
sin 2 0 23 


This dependence on 0 looks complex, but the boundaries 
of the 5 = 0 regions can be seen in Fig. [3] to be those val¬ 
ues of 0 for which one mixing angle is either 0° or 90°, and 
closed form expressions can be derived for these values 
from the mixing angle formulae above. For the matrix 
shown here, this means that only the S = ir solution is 
agrees with the global data. However, when considering 
all permutations the alternative CP conserving solution 
can also be found. 

The prediction for 0 23 can be expressed as an atmo¬ 
spheric sum rule to first order in r, 


l 2 _ 1 + -| + 2 r + °( r2 )' 

1 + v 5 1 T 


This relation, being derived from a non-cyclic symme¬ 
try in the charged-lepton sector, has to the best of our 
knowledge not been presented before in general analyses 
of atmospheric sum rules m- 

A permutation in the 023 plane acting from the left of 
the PMNS matrix effects a mapping of 023 —> f — 023- For 
this reason, both of the expressions defined above have 
a complementary pattern with an inverted 023- These 
alternative patterns are shown by dashed lines in Fig. [3] 


F. Predictions with two and three degrees of 
freedom 

So far we have analysed the cases when G e £ 
{Z 3 ,Z 5 ,Z 2 x Z 2 } and G v = Z 2 . The patterns resulting 
from these groups have a single degree of freedom, the 
unphysical angle 0 , controlling their mixing parameter 
predictions. There are however, more general scenarios 
where the 3 angles and 3 phases of the PMNS matrix are 
specified by 2 or 3 input parameters. 

The cases with 2 degrees of freedom arise from the 
choice G e = Z 2 while the neutrino symmetry is enlarged 
to the full Klein group, G„ = Z 2 x Z 2 . The symmetry of 
the charged-lepton mass terms is insufficient to uniquely 
specify the diagonalising matrix of the mass matrix, and 



Internal 9 (deg.) 

FIG. 3. Allowed mixing angles for G e = Z 2 x Z 2 as a function 
of the unphysical parameter 9. There are two possible sets of 
predictions of the mixing angles which have the same 0 i 2 and 
0 i 3 predictions but distinct 023 predictions (solid and dotted 
lines) related by the mapping 623 f — 023- All complex 
phases are CP conserving for these patterns: sin 5 = sin Q 21 = 
sin 031 = 0. The shaded regions show the 3<r allowed region 
for the corresponding mixing angle according to current global 
data |25| . 


requires in general a further 2 -dimensional complex rota¬ 
tion, 


I 1 0 0 . \ 

i? e ( 0 ,y) = 0 cos 0 sin 0 e 17 I . 

\0 — sin 0 e -17 cos 0 J 

However, the neutrino symmetry has been enlarged, and 
so t/„ is uniquely specified by the symmetry generators. 
We have scanned over all such combinations and found 
that no viable patterns arise from this scenario. 

There is one more combination of residual symmetries 
possible in our construction: G e = Z 2 and G„ = Z 2 . This 
is an extension of the previous case, where the charged- 
lepton symmetry introduces two parameters but now the 
neutrino residual symmetry also requires a single real pa¬ 
rameter to diagonalise the most general mass matrix. In 
principle, there is no reason to discount these patterns. 
They are consistent with the idea that the full flavour 
group has broken into residual subgroups implying cor¬ 
relations on the flavour observables. However, the in¬ 
creased number of parameters reduces the predictivity of 
the theory (3 inputs, 6 outputs). We have scanned over 
such groups and verified that there are eligible patterns 
which match all the global data. Some of these patterns 
feature non-constant Majorana and Dirac phase predic¬ 
tions, and many are not simply related to the patterns 
that we have found in more restrictive schemes. However, 
due to the larger parameter space and reduced predictiv¬ 
ity, we will not attempt to present any results of this 
type. 
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IV. PHENOMENOLOGICAL PROSPECTS 

In the preceding sections we have derived all patterns 
of mixing angles and phases which are possible with an 
A 5 symmetry with generalised CP broken into residual 
symmetries. They depend upon a single real angle, 0, 
and can all be brought into agreement with current global 
data |25] for a suitable restriction of its range. Eliminat¬ 
ing the unphysical parameter 8 leads to a set of corre¬ 
lated predictions between observables which are testable 
by oscillation experiments and searches for neutrinoless 
double-beta decay. In this section, we shall discuss the 
prospects for present and future experiments to constrain 
these patterns and derive simple versions of the predicted 
parameter correlations which may be useful experimen¬ 
tally. We would like to stress that the correlations iden¬ 
tified in this paper will be tested at almost every stage in 
the experimental programme of the next few decades. 
Near term results from T2K m and NO^A [58! can 
be expected on the maximality of 9 23 and S, and in the 
medium term, new reactor and long-baseline experiments 
such as JUNO [59], RENO-50 [60], DUNE^ T2HK [H] 
and possibly ESS^B [62] should bring us increased pre¬ 
cision on 812 , 823 and S. Finally, for the most stringent 
tests of the models in question, the option remains to 
construct a more ambitious facility such as the Neutrino 
Factory [65] . In a complementary direction, neutrinoless 
double beta decay experiments will further sensitivity to 
this decay, providing evidence on the Majorana nature 
of neutrinos and, at least in principle, the first measure¬ 
ments of the values of Majorana phases. Many of these 
observations are largely independent and we can expect 
significant evidence either in favour of, or ruling out, the 
patterns identified in this paper. 


A. Precision measurements of 6 12 

The viable sets of mixing parameters which we have 
found above predict correlations in 9 12 and 6 * 13 , and 
therefore very precise measurements of these angles have 
the potential to discriminate between flavour symmetric 
patterns, or to rule them out entirely [23 . 

The upcoming medium-baseline reactor (MR) neutrino 
oscillation experiments, such as JUNO 155] and RENO- 
50 [60], expect to make very precise, sub-percent mea¬ 
surements of the oscillation parameter 0 12 . The preci¬ 
sion on 0 i 3 , currently dominated by measurements from 
Daya Bay [2] and RENO [2], is not expected to be signif¬ 
icantly improved by the next generation of reactor facili¬ 
ties. Therefore, the first significant test of the predictions 
of this paper will come from increased precision on 0 i 2 
independently of (A 3 . We have identified 3 distinct pre¬ 
dictions for 0 i 2 , if we fix 0i 3 to its current best-fit [25] 


5 The new name for the LBNF/ELBNF project. 



6,3 (deg.) 


FIG. 4. Predictions for 0i2 as a function of 0i3. All of 
the patterns of mixing parameters associated with a given 
charged-lepton residual symmetry have the same prediction 
(solid lines). The dashed line close to each prediction shows 
the linearised predictions in Eq. (151, Eq. (161 and Eq. (171. 
The grey regions show the 1 and 3cr allowed regions for 0i2 
from current global data [25 ;. 


these are 


012 = 35.71°, 012 = 32.11°, 0 12 = 35.14°, 


for preserved charged-lepton subgroups Z 3 , Z 5 and Z 2 x 
Z 2 , respectively. Given that the expected precision of the 
MR experiments is at the level of 0.1° or around 0.3% 
for 0 i 2 , a strong discriminatory power exists between the 
values of the mixing angles predicted by these correla¬ 
tions. The difference between the predicted values of all 
models under consideration is always greater than 0.26° 
over the current 3er interval for 0 i 3 , and in many cases 
significantly greater. Therefore we can expect these ex¬ 
periments to identify with considerable confidence if any 
of the charged-lepton residual symmetries are consistent 
with observation. 

In the framework discussed in this article, each model 
predicts a continuous correlation between the values of 
013 and 0 i2 . If one of the predictions above appears to 
agree with data, it would be desirable to test the correla¬ 
tion between parameters itself. These correlations can be 
conveniently expressed as expansions in the dimension¬ 
less parameter r = v / 2 sin 0 i 3 m- The current global 
best-fits give 0 i 3 « 8.50° [25] which translates to r ~ 0 . 2 ; 
the second-order corrections are therefore suppressed by 
a factor of 1/25. Expressed in this way, the predictions 
for sin 0 i2 associated with the charged-lepton subgroups 
Z 3 , Z 5 and Z 2 x Z 2 (respectively) can be expanded in the 
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following relations, 



Expressing these in terms of the angles themselves, we 
find 


012 = 35.27° + 10.13° r 2 + 0 (r 4 ) , (15) 

0i 2 = 31.72° + 8.85° r 2 + O (r 4 ) , (16) 

0i2 = 36.00° - 19.72° r 2 + O (r 4 ). (17) 


The approximations in Eqs. (15), (16) and have been 
plotted against the unapproximated expressions for 0 12 
in Fig.|4j We see that these relations depend only slightly 
on 0 i 3 , which first appears at the order (D(r 2 ), leading to 
sub-degree level corrections. 

The formulae above show that the predictions for 0 12 
only vary by 0.07°, 0.06° and 0.13° (for Z 3 , Z 5 and 
Z 2 xZ 2 , respectively) over the current 3cr region for 0i 3 . 
This is of the order of the target precision of the MR 
experiments, and it is therefore unlikely that the 0 12 - 
0i 3 correlations themselves will be tested at a significant 
level even if precision on 0i 3 were to be greatly improved. 
There are no currently planned facilities which could fur¬ 
ther improve the precision on 0 i 2 - 


B. Maximal-maximal predictions for 623 and S 


The current and upcoming generation of long-baseline 
experiments will be able to place important constraints 
on the parameters 02 3 and S. Measuring 02 3 and <5 inde¬ 
pendently will provide valuable information on the via¬ 
bility of flavour symmetric models; however, in the pat¬ 
terns that we have identified the maximality of 02 3 is 
significantly correlated with the value of 5. In four of 
these patterns (excluding different Majorana phase as¬ 
signments), two from G e = Z 3 and two from G e = Z 5 , 
predict a maximal value of 02 3 and a maximal amount of 
CP violation, 


023 — 


7r 
4 


and 


|sin<5| = 1. 


Therefore, the joint determination of these parameters 
around these maximal values would be a particularly in¬ 
teresting measurement from the point of view of GCP 
model building. 

Testing the maximality of these parameters is an at¬ 
tainable goal for current and future oscillation experi¬ 
ments. After its full period of data taking, T2K expects 
to be able to exclude maximal 02 3 at the 90% C.L. for 
Isin 2 (20 23 ) — 0.5| > 0.05-0.07 largely independently of 


the value of 5 m- Measuring S itself is significantly 
harder; however, the maximal CP violating values con¬ 
sidered here are the most accessible. T2K can expect to 
be able to exclude 0<<5<7 t(7t<6< 27t) at the 90% 
C.L. for a true value of d = 3tt/2 (<5 = tt/ 2 ) [FT;. This 
would allow T2K to distinguish between 5 = n/2 and 
5 = 37 t/2 if one of them is true at at least the 90% C.L. 
NOi^A can also be expected to contribute to this mea¬ 
surement |58j with a similar power for excluding S = tt/2 
and 5 = 37 t/2 . Although these exclusions are expected 
to be individually statistically weak, they would consti¬ 
tute valuable information on the validity of the models 
studied in this paper and would provide strong encour¬ 
agement for future work by the next-generation of oscil¬ 
lation experiments. 

In the medium term, new long-baseline experiments 
are expected with significantly improved sensitivities, in 
particular to the phase S, allowing the maximal-maximal 
predictions to be further tested. To estimate the po¬ 
tential for excluding these models with future facilities, 
we have run a simulation of LBNE/DUNE using the 
GLoBES package |(M|. Our simulation is based on the 
detector responses files and fluxes made available by the 
LBNE collaboration in Ref. US)- We point out that 
thanks to its more ambitious design, it seems likely that 
the DUNE project can significantly improve the sensitiv¬ 
ity computed here. However, without access to updated 
experimental information, making a quantitative assess¬ 
ment of the extent of this improvement is challenging. We 
assume a 700 kW beam operating at 120 GeV, a detec¬ 
tor based on liquid Argon-TPC technology with a mass 
of 34 kton, and overall systematic errors of 5% for both 
the signal and background normalizations. The results of 
these simulations are shown in Fig. [5j where we present 
the regions of true parameter space for which the com¬ 
binations of (023 ,5) = (7r/4, 7 t/2 ) and (7r/4,37r/2) can be 
excluded after 5 years neutrino and 5 years antineutrino 
running. We find that these patterns can be excluded 
at 3 a if the true value of 0 23 satisfies 9 23 < 43.0° or 
023 ^ 48.3°, or if the true value of S is outside the in¬ 
tervals 90°+6go or 270°1 ^o. Here we see the importance 
of observing both 023 and S for excluding our models. 
A measurement of 9 2 3 alone will not be able to distin¬ 
guish between the models which predict 023 -maximality 
and the model from Z 2 x Z 2 which predicts values of 02 3 
which differ from maximality by only around 2°. How¬ 
ever, these models have maximally distinct predictions 
for <5, and as we have shown, the measurement of <5 alone 
would be able to separate these cases at 3 ct. 


C. Dirac CP conserving patterns and precision 
measurements of 023 

Those patterns which instead make non-maximal pre¬ 
dictions of 023 also predict CP conserving values of <5, 
such that | cos S | = 1, and we can expect constraints 
to be placed on these models by the attempts to dis- 
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FIG. 5. Red (green) lines show the exclusion regions at 1, 2 
and 3cr for #23 = 7r/4 and <5 = 37 t/ 2 (<5 = 7r/2) expected at 
LBNF with a 34 kton LAr detector after 5 + 5 years running. 
In this regions outside the curves, the two sets of predictions 
can be excluded at the given confidence. The side panels 
show the appropriate marginalised Ay 2 and the 1, 2 and 3<r 
confidence levels (1 d.o.f). 


cover leptonic CP violation — a standard search for the 
next-generation of CP-sensitive oscillation experiments 
[61] l! 62l f66) - l69| . It has been shown that LBNO running 
with a beam derived from the SPS accelerator at CERN 
could rule out leptonic CP conservation at 3er for around 
45% (65%) of the parameter space for a detector mass of 
20 kton (70 kton). This could be increased to 70% (80%) 
with an upgraded beam power m- LBNE has predicted 
a similar sensitivity to CP violation |69j . with the ulti¬ 
mate reach also depending crucially on the planned series 
of upgrades to detector mass and beam power. With a 
10 kton detector and 6 years of data using a 1.2 MW 
beam, the measurement could be made for 33% of the 
parameter space at 3<r. This rises to 40% of the parame¬ 
ter space at 5er once the detector mass has been increased 
to 34 kton and 6 years more data has been collected. Fi¬ 
nally, a beam power upgrade to 2.3 MW could increase 
this to 60% of the parameter space at 5a [69]. The T2HK 
and ESS^B proposals also show strong sensitivity to CP 
violation, both using a megaton-scale water Cerenkov de¬ 
tector and MW power beams. T2HK has shown that it 
can expect a discovery of CP violation over 76% (58%) 
of the parameter space at 3er (5cx) |70J • A similar reach is 
possible with ESS^B, which expects a 3er (5a) discovery 
of CP violation after 10 years of data-taking over 74% 
(50%) of the parameter space [62 1. 

If the current experimental programme fails to discover 
CP violation in the leptonic sector and sin S is discovered 
to be small, there are four distinct patterns from our 
model which would remain in agreement with the data. 
These models can be tested by the increased precision 
on measurements of #23 expected from next generation 
long-baseline facilities. In all cases of this kind, we find 
predictions coming in pairs. The model associated with 



FIG. 6. #23 as a function of #13 for the patterns which predict 
CP conservation (solid lines). The dashed line close to each 
solid line shows the linearised expression in Eqs.( |18| l and ( | 19| ). 
The grey regions show the 1 and 3cr allowed regions from 
current global data [25] . 


a Z 5 residual symmetry predicts 

0 23 = 45° ± 25.04° r + 0(r 2 ), 


(18) 


while the predictions for a residual Z 2 x Z 2 symmetry are 
given by 


# 23 = 31.72° + 55.76° r + 0(r 2 ), 
#23 = 58.28° - 55.76° r + 0(r 2 ). 


(19) 


The pairs of sum rules given by Eq. (18) or Eq. (19) 


are related by the octant degeneracy, #23 —> § — #23, as 
can be seen clearly in Fig. [ 6 ] where the approximations 
above are shown against the full predictions. We note 
that in contrast to those for # 12 , these relations depend 
on r at linear order. Therefore, they are far more sensi¬ 
tive to the precise correlation between parameters, and 
the measurement of the correlation itself becomes more 
accessible. 

The first discriminating factor between these solutions 
will come from improved precision on #23 and the resolu¬ 
tion of the octant degeneracy. In these models we predict 
non-maximal mixing, and a successful determination of 
the octant would provide early evidence in their favour. 
This would be most challenging for the model based on 
Z2 x Z2 for which #23 differs from 45° by between 2.5° and 
0.8° over the current 3cr range of #13. The model based 
on Z 5 instead predicts greater deviations from maximal 
atmospheric mixing, ranging between 4.8° and 5.6° over 
the same interval. Studies of the potential for the current 
generation of oscillation experiments, of which T2K and 
NOi'A play the most important role, suggest that the 
octant can be established at 3<r (2er) for deviations from 
maximality greater than around 6 ° (4°) [71] [Hi. This 
precludes the current generation from separating between 
the two predictions of Z2 x Z2, but would allow for 2 a evi¬ 
dence for those predictions coming from our model based 
on Z5. This discovery potential will be improved by the 
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next generation of oscillation experiments. In Ref. ]69j it 
is shown that with an exposure of 60 kton-years, LBNE 
could determine the octant at 3cr if the true value of $23 
deviates from maximality by more than 4°-5°. However, 
the best bounds could come from T2HK by studying at¬ 
mospheric neutrino data. A 3er determination of the oc¬ 
tant is expected to be possible after 10 years of data- 
taking for true values |sin 2 023 — 0.51 > 0.04-0.06 cor¬ 
responding to deviations between 2°-3° m- Although 
exclusion of the Z 2 x Z 2 pattern would be unlikely, the 
two predictions from Z 5 would be distinguishable. 

To go beyond the octant measurement, higher preci¬ 
sion will be necessary to separate between the Z 5 and 
Z 2 x Z 2 predictions, or indeed to test their specific corre¬ 
lations with 0i3. The difference between these two pre¬ 
dictions varies from between 2.4° to 4.8° over the current 
allowed regions. Therefore degree-level precision will be 
required to distinguish between them, even in the pres¬ 
ence of greatly improved knowledge of 0i3. In Ref. [69] it 
is shown that the minimal 10 kton LBNE configuration 
running for 6 years would have a precision of around 1° 
at la for true values of 023 around 51°, which increases 
as we approach 023 maximality to a lcr width of around 
2.5°. Similarly, T2HK shows that around the point ex¬ 
pected to give the worst sensitivity to 023, the 90% C.L. 
width is around 2°-3° m- These results suggest that 
a significant discrimination between these models would 
be challenging with these set-ups; however, evidence in 
favour of these models would be possible at low signifi¬ 
cance and if observed in conjunction with an absence of 
observable CP violation, this would present a concrete 
hypothesis for future work. 


D. Long-term prospects 

We have seen in the previous sections that although the 
next generation of superbeam and reactor experiments 
will be able to test the consistency of the patterns that 
we have identified, much of their discriminatory power 
relies on excluding maximal angles and phases. Testing 
the continuous correlations predicted in our models, for 
example between 0 12 and 0 i 3 or between 0 2 3 and 0 i 3 , 
would require higher precision. 

The only proposed experiment capable of pushing 
the precision frontier beyond the results of the next- 
generation superbeams is the Neutrino Factory (NF) [73] . 
which produces a beam with low systematic uncertain¬ 
ties from the decay of stored muons [55] . A NF would 
be able to improve our knowledge of the mixing pa¬ 
rameters in a number of ways, but for the present pur¬ 
poses it serves two main roles. Firstly, such a facility 
would greatly increase the precision on d, with an ulti¬ 
mate lcr precision estimated at around 5° [7U [75]. This 
could allow many of our models to be excluded indepen¬ 
dently of their other parameter correlations. Secondly, 
a NF would provide a high-precision determination of 
$23 m, allowing it to perform a very stringent test of 



FIG. 7. The regions of true parameter space for which the 


upper-octant relations in Eqs. (181 and (191 can be excluded 
at 1, 2 and 3 ct by a 2000 km, 10 GeV LENF using a MIND 
detector. The two regions start to overlap at low values of 0i3 
at the 5 ct confidence level. 


the 023-013 correlations discussed in the previous subsec¬ 
tion. To quantify this possibility, we have performed a 
simulation of a representative Low-Energy Neutrino Fac¬ 
tory [73 [77] to compute the regions of true parameter 
space which would allow two candidate models to be (in¬ 
dividually and collectively) excluded. Our simulations 
assume a 2000 km baseline and a stored-muon energy 
of 10 GeV which is close to the optimal configuration 
for CP violation discovery [75H5T| . For our detector, 
we take a Magnetized Iron Neutrino Detector (MIND) 
[82l [83] , assumed to be in a toroidal magnetic field allow¬ 
ing for muon charge identification, with a fiducial mass of 
100 kton. The detector response is described by a set of 
migration matrices provided by Ref. [84j . Backgrounds 
to the appearance channel signal come from both charge- 
and flavour-misidentified events, as well as the secondary 
decay products arising from t± decays in the detector. 
An overall 1% (10%) uncertainty is taken on the signal 
(background). Our results are shown in Fig. [ 7 } where the 
blue (red) shaded regions show the area of true param¬ 
eter space for which the Z 2 x Z 2 (Z 5 ) correlation would 
not be able to be excluded. We see that the two coloured 
regions do not overlap, and therefore all points in this 
parameter space allow for the exclusion of at least one of 
the correlations at 3cr or higher significance. 


E. Neutrinoless double beta decay 

Over the next decade, the new generation of neutri¬ 
noless double beta (Oi//3/3) decay experiments will signif¬ 
icantly increase the sensitivity to this rare process. For 
the first time these experiments will probe the region of 
parameter space associated with the inverse hierarchi¬ 
cal spectrum. These experiments aim to establish that 
neutrinos are Majorana in nature, but can also provide 
valuable information on the neutrino mass spectrum and 
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in principle, measure the Majorana phases themselves. 

The CK/3/3 decay rate is proportional to the effective 
Majorana mass \m ee \ (see e.g. Ref. J5TJ E5|), which is 
given by 


m Ke \ = 


U ek m k 


1 


mi cos 2 0 12 cos 2 $13 + m 2 sin 2 $12 cos 2 Oise' 1 ' 


+m 3 sin 2 0 i 3 e ^“ 31 2<5) 


( 20 ) 


where ol^i and 0(31 are Majorana phases and 6 is the Dirac 
phase. The predicted values of \m ee \ depend crucially on 
the neutrino masses. The latter can be ordered in two 
ways: normal ordering (NO; mi < m 2 < m 3 ) or inverted 
ordering (10; m 3 < m\ < m 2 ). As the parameters Am^ 
and Am 2 -, are known from oscillation physics, there is a 
single degree of freedom remaining amongst the masses. 
This is typically taken to be the lightest neutrino mass, 
mi (m 3 ) for NO (10), which we will denote in both cases 
by mi [85]. The parameter space available to |m ee | can be 
further divided into three particularly interesting regions 
based on the true value of mi. The first is for quasi¬ 
degenerate masses (QD) where mi > 0.1 eV, in which 
the splitting between masses is a small correction to ap¬ 
proximately degenerate values. For smaller values of mi 
there are two parameter regions: one for normal hierar¬ 
chical masses (NH; m\ < m 2 m 3 ) and the other for 
inverted hierarchical masses (IH; m 3 mi < m 2 ). In 
order to better understand the predictions of our models, 
we have first computed the predicted values of |m ee | in 
the generic case, assuming only that the mixing param¬ 
eters lie in their current 3cr allowed ranges [25]. These 
predictions for NO (10) are shown as the blue (red) re¬ 
gion in Fig. [ 8 ] For quasi-degenerate and IH spectra, there 
exist lower bounds on |m ee | [ 86 ]: for IH |m ee | > 0.015 
eV, and for QD |m ee | > 0.03-0.04 eV. For NO there is 
no non-zero lower bound as \m ee \ can vanish due to a 


cancellation between terms in Eq. (20). 


There are many experiments that are searching for 
On/3/3 decay or are in various stages of planning and 
construction. The most recent experiments that have 
set upper limits on the effective Majorana mass are 
CUORICINO, GERDA, EXO-200 and Ea r n LAND-Zen. 
CUORICINO, based in Gran Sasso National Laborato¬ 
ries in Italy, was an experiment used to test the feasi¬ 
bility of its successor, CUORE. Using data taken from 
2003-2004 they achieved a bound on |m ee | of 200-1100 
meV EZ1, where as with all Oz^/3/3 experiments a key un¬ 
certainty on their limits comes from the nuclear matrix 
element. GERDA (The Germanium Detection Array) 
is also located in the Gran Sasso Laboratory LNGS in 
Italy [55] ■ During phase I of their data taking period 
they acquired sufficient data to attain an upper limit of 
|m ee | < 200-400 meV. They intend to increase their sen¬ 
sitivity by a factor of approximately ten during GERDA 
Phase II. EXO (Enriched Xenon Observatory)-200, lo¬ 
cated in Carlsbad, New Mexico, has placed upper bounds 


for |m ee | of 69-163 meV [59]. KamLAND-Zen (Kamioka 
Liquid Scintillator Anti-Neutrino Detector-Zero neutrino 
double beta decay search) has had two phases of data 
acquisition and using the combined data they found an 
upper limit of |m ee | < 140-280 meV [99] ■ Both EXO- 
200 and KamLAND-Zen have ambitious long term plans 
to upgrade their experiments in order to explore the 
inverted hierarchical region of parameter space. EXO- 
200 intends to upgrade to nEXO (next Enriched Xenon 
Observatory) which, with ten years of data taking, is 
expected to cover the full IO region m- In its next 
phase, KamLAND-Zen aims to increase its sensitivity to 
around 50 meV after approximately two years running 
time (90). As this limit will only begin to probe the IH 
region, KamLAND-Zen has proposed KamLAND2-Zen. 
This upgraded detector (with a running time of 5 years) 
has a target sensitivity of |TO ee | — 20 meV and will allow 
the exploration of the majority of the IO region and all of 
the QD parameter space. In addition to EXO-200 and 
KamLAND-Zen, other future 0^/3/3 decay experiments 
are CUORE, SNO+ and NEXT. CUORE plans to start 
taking data in 2015. Over the course of five years, they 
hope to reach a sensitivity of 50-120 meV [92] , SNO+, 
a multi-purpose experiment located in Sudbury, Canada, 
aim to achieve a similar upper bound on |m ee |. After 
two years of data taking they expect to be able to set the 
upper bound |m ee | < 100 meV [93]. The NEXT (Neu¬ 
trino Experiment with Xenon TPC) experiment, based 
at CanFranc Underground laboratory (LSC), will com¬ 
mence data taking in 2018 using the NEXT-100 detector. 
Despite their late start compared with that of EXO-200 
and KamLAND-Zen, they intend to achieve a |m ee | sen¬ 
sitivity of approximately 100 meV by 2020 [94]. The next 
stage will be the development of BEXT which proposes 
to fully cover the predicted values of |m ee | for IO [95] , 
Although they are not discussed here, there are other ex¬ 
periments which aim to improve the current bounds on 
|m ee |; for example, COBRA [99] . the Majorana Demon¬ 
strator [97] , SuperNEMO [95] and the DCBA experiment 
99.. amongst others. 

The predicted values of |m ee | for the case of A 5 with 
GCP can be calculated from the leptonic mixing matri¬ 
ces of Eq. (10), Eq. (11) and Eq. (14) for both IO and 
NO. The complex phases only influence |m ee | through the 
combinations e la21 and e^ a31 ~ 2S \ and we will denote the 
phases of our predictions by an ordered pair of ± signs 
e.g. (- 1 —) when 0:21 = 0 and « 2 i — 26 = it. As |m ee | does 
not depend on 6 * 33 , patterns which only differ by this an¬ 
gle will be degenerate and each preserved charged-lepton 
subgroup leads to a single prediction for each mass or¬ 
dering and phase assignment. Fig. [ 8 ] shows the predicted 
values from the mixing patterns in this paper for each 
charged-lepton residual symmetry G e . In these plots, we 
have neglected a small width to each line which comes 
from varying 0 13 and the neutrino mass-squared split¬ 
tings over their allowed ranges, instead fixing these at 
their best-fit values from Ref. ) 25] , 

We focus first on the IO spectra. The phase assign- 
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FIG. 8 . \m ee \ versus the lightest neutrino mass for the 10 (NO) for the solid (dashed) lines. The predictions in a given panel 
all have the same phase assignment, shown in the top left of the plot. The red (blue) shaded region shows the most general 
predictions for |m ee | with IO (NO) obtained by varying the oscillation parameters over their current 3cr global ranges f25l j. 


ments (++) and (d—), shown on the top row of Fig. | 8 j 
predict large values of \m ee \, close to the upper bound¬ 
ary of the 10 region obtained using the 3 a global data. 
These predictions are very similar for all models. This 
can be understood as the term in \m ee \ proportional to 
m 3 only has a subdominant effect: it is not only mul¬ 
tiplied by the small number sin 2 0 13 , but is further sup¬ 
pressed for IH by the small value of m 3 itself. If we 
neglect this term, the resulting approximation at leading- 
order is independent of 0 i 2 up to corrections of the order 
OfAniji/Amjj). It is feasible that experiments such 
as CUORE and KamLAND-Zen, and to a much greater 
extent BEXT and nEXO will be able to explore this top¬ 
most region of the IH parameter space and test these 
predictions. Further distinguishing between them will 
be beyond their scope due to the small predicted differ¬ 
ences and substantial experimental and theoretical un¬ 
certainties on \m ee \. For the phase assignments (—b) 

and (-), shown on the bottom row of Fig. [ 8 j we see 

values of \m ee \ that are further suppressed and which 
exhibit more model dependence. Once again, the sup¬ 
pression of the m 3 term explains the similarity between 
the two phase assignments. The lower values compared 
to the (++) case arise from the relative phase differ¬ 
ence between the m\ and m 2 terms: at leading order 


|m ee | = -^/fArnlTf cos 2 0i 3 cos 2 ( 2012 ) [5S11100] . This ef¬ 
fect is evident if we compare Z 5 and Z 3 cases: the larger 
012 value of Z 3 accounts for the more pronounced cancel¬ 
lation and therefore lower \m ee \ than that of Z 5 . These 
predictions are beyond the reach of many of the facilities 
discussed so far; although KamLAND2-Zen, aims to set 
limits near the predictions for Z 5 and, if capable of test¬ 
ing the full IO region, should be accessible to nEXO and 
BEXT. Although lying in a region of parameter space 
that is harder to explore, the greater model dependence 
for these phase assignments would make it easier to dis¬ 
tinguish between models than with the (++) and (H—) 
cases. There exists a separation of around 5 meV be¬ 
tween the predictions for Z 5 and the other subgroups; 
however, it is unlikely such a resolution on \m ee \ would 
be attainable in the foreseeable future. 

For NO, we see quite different behaviour. In the quasi¬ 
degenerate region, the mass-squared splittings are negli¬ 
gible and the predictions for the IO and NO cases effec¬ 
tively coincide. However, in the limit of vanishing m; the 
situation is very different. In this limit it is the relative 
phase between the m 2 and m 3 terms which dominates 
the magnitude of |m ee |, which leads to larger predictions 
for the phase assignments (++) and (-), while sup¬ 

pressing the predictions of (—b) and (H—). Although 
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G e 

012 

023 

sin a.ji 

5 

Z 3 

35.27° + 10.13° r 2 

45° 

0 

90° 

270° 



45° ± 25.04° r 

0 

O 

O 

Z 5 

31.72° +8.85° r 2 

180° 

45° 

0 

90° 



270° 



31.72° + 55.76° r 

0 

O 

O 

Z 2 X Z 2 

36.00° - 34.78° r 2 

180° 


58.28° - 55.76° r 

0 

O 

O 



180° 


TABLE I. Numerical predictions for the correlations found 
in this paper. The dimensionless parameter r = y/2 sin #13 is 
constrained by global data to lie in the interval 0.19 < r < 
0.22 at 3 a. The predictions for 6*12 and 023 shown here ne¬ 
glect terms of order 0(r 4 ) and 0(r 2 ), respectively. Following 
the method of this paper, the Majorana phases can only be 
predicted modulo 7 r and the values in the fourth column hold 
for all phases. 

exploring the NH region experimentally is beyond the 
scope of any planned experiment, if 0i//3/3 decays are not 
observed and oscillation physics establishes that the neu¬ 
trino masses are NO, it would be of paramount impor¬ 
tance to try and test \m ee \ values in the NH region. Due 
to the rich interplay between relative phases, these mod¬ 
els make quite different predictions across this parameter 
space. In fact, all mixing angle patterns discussed in this 
paper could accommodate a value of \m ee \ near the top of 
the current NH region allowed by global data. Although 
such an observation would add further support to any 
prediction of this paper which was still consistent with 
experimental data, to further discriminate between these 
models it would be necessary to provide complementary 
information on the absolute mass scale. 


V. CONCLUSIONS 

Assessing the viability of flavour symmetric models of 
the leptonic sector is an accessible target for precision 
measurements from present and future neutrino oscilla¬ 
tion experiments. In this article, we have presented a 
detailed analysis of a particular theoretical scenario: the 
flavour symmetry A 5 with a generalised CP symmetry 
breaking into residual subgroups at low energies. We 
have identified the most general form of the generalised 
CP transformation, and studied the full group for con¬ 
sistent residual symmetries. Our analysis results in 6 
distinct sets of mixing angle predictions each with an 
additional 8 possible combinations of phases which are 
shown in Table [Tj These depend at most on a single real 
parameter, and predict testable correlations between cer¬ 
tain parameters. In addition, the Majorana phases for all 
of our predictions are CP conserving. These patterns can 
be classified by the residual symmetry in the charged- 


lepton mass terms: Z 3 , Z 5 and Z 2 x Z 2 . A symmetry 
of Z 3 predicts maximal 023, maximal CP violation from 
S and a value of 0 i2 that lies close to the upper bound¬ 
ary of the 3 a global fit data. There are two distinct pat¬ 
terns which arise from a preserved Z 5 residual symmetry. 
These share a common $12 prediction which lies close to 
the lower boundary of the 3tr global fit data; however, 
one prediction has maximal 0 23 and a maximally CP vi¬ 
olating value of 5 whilst the other has non-maximal 023 
and CP conserving values of <5. The patterns arising from 
a preserved subgroup Z 2 x Z 2 also share a common 0 12 
which lies above the current la region. In this case both 
023 predictions are non-maximal and the value of S is 
CP-conserving. 

We have then discussed the phenomenology of our pre¬ 
dictions, focusing on the role which current and future 
reactor, superbeam and neutrinoless double beta decay 
experiments can play. The predictions for 0 i 2 should be 
testable at high significance by the next generation of 
reactor neutrino experiments, such a JUNO and RENO- 
50. These experiments can be expected to distinguish 
between the different models; however, testing the pre¬ 
cise correlations between 0 12 and 0 13 will most proba¬ 
bly remain beyond the reach of any foreseen experiment. 
A particularly interesting feature of the patterns found 
in this paper is the correlated maximality of 023 and <5, 
and also non-maximal 023 and CP conserving values of 
5. Testing these correlations is a feasible goal for current 
and future superbeam experiments. T2K and NOi^A can 
be expected to collect early evidence if such a pattern ob¬ 
tains, and we have shown that DUNE will be able to iden¬ 
tify such a pattern over a significant part of the parame¬ 
ter space. For the CP conserving patterns, the deviations 
from 023 = 7 r /4 are expected to be measureable at 3<r by 
the next generation of superbeams for the preserved sub¬ 
group Z5, but not for Z2 x Z2. Ultimately separating 
between these models at 3er significance across the whole 
parameter space could be done using a Neutrino Factory 
after 10 years of data taking. An attractive feature of the 
theoretical scenario in this work is its ability to predict 
Majorana phases, and therefore, observables for neutri¬ 
noless double beta decay experiments. We have seen that 
in the case of inverted mass ordering, two of the possi¬ 
ble Majorana phase combinations predict the discovery 
of neutrinoless double beta decay at upcoming experi¬ 
ments. In the longer term, the exploration of the full 
parameter space for inverted hierarchical mass spectra 
could allow all of our patterns with this mass spectrum 
to be confirmed independently of oscillation physics. 

In conclusion, we find that the combination of the 
flavour symmetry A 5 with a generalised CP symmetry 
allows for a number of viable predictions to be made for 
the mixing angles and phases. These predictions spec¬ 
ify parameter correlations which present good targets for 
each stage of the next decade of the experimental pro¬ 
gramme. 

Note added: During the final preparations of this ar¬ 
ticle, preprints of two similar works were made available 
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mm- These works also study the mixing patterns 
arising from the residual symmetries of A 5 with GCP. 
The patterns derived in the first part of our paper con¬ 
firm those found in Refs. 11011 and [1021 . Although, case 
II in Ref. [lOlj is omitted in our analysis as its predic¬ 
tions fall outside the 3cr global intervals used in this work. 
Our phenomenological work, however, significantly ex¬ 
tends the analysis in these papers. 
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